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Kondo insulators are particularly simple type of heavy electron material, where a filled band of heavy quasi- 
particles gives rise to a narrow band insulator. Starting with the Anderson lattice Hamiltonian, we develop 
a topological classification of emergent band structures for Kondo insulators and show that these materials 
may host three-dimensional topological insulating phases. We propose a general and practical prescription of 
calculating the Z2 strong and weak topological insulator indices for various lattice structures. Experimental 
implications of the topological Kondo insulating behavior are discussed. 
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Kondo insulators are a particularly simple type of heavy 
fermion material, first discovered forty years ago[ 1], in which 
highly renormalized f-electrons, hybridized with conduction 
electrons, form a completely filled band of quasiparticles with 
excitation gaps in the millivolt range ylyQ. While these ma- 
terials are strongly interacting electron systems, their excita- 
tions and their ground-states can be regarded as adiabatically 
connected to non-interacting band-insulators @|. 

It was recently shown that time-reversal invariant band in- 
sulators can be classified by the topological structure of their 
ground-state wavefunctions |S la [l |8[] . One of the dramatic 
consequences of this discovery, is the existence of a new 
class of "topological" band-insulator in which strong spin- 
orbit coupling leads to a ground-state that is topologically dis- 
tinct from the vacuum, giving rise to gapless surface excita- 
tions. 

In this paper, we show that Kondo insulators, as adiabatic 
descendents of band insulators, can also be topologically clas- 
sified. The strong spin-orbit coupling characteristic of these 
materials leads us to predict that a subset of Kondo insulators 
are topologically non-trivial, with anomalous surface excita- 
tions. In current models of topological insulators, the spin- 
orbit coupling is encoded in a spin-dependent hopping am- 
plitudes between different unit cells. By contrast, in a topo- 
logical Kondo insulator (TKI), we show that the topologically 
nontrivial insulating state is produced by the spin-orbit cou- 
pling associated with the hybridization between conduction 
and /-electrons. 

Below, we develop a simple model for topological KIs. The 
physics we study is motivated by the canonical Kondo insulat- 
ing behavior of SmBell] andCe3Bi4Pt3 [9]. The realization of 
a particular topologically nontrivial insulating state depends 
on the position of renormalized /-level relative to the bottom 
of the conduction band (see Fig. 1). To analyze the topology 
of the bands in these materials, we use a periodic Anderson 
lattice model. 

In a KI, the insulating state arises due to hybridization be- 
tween the conduction and /-electrons, provided that the chem- 
ical potential lies inside the hybridization gap separating the 
quasiparticle bands. The spatial symmetry of the hybridiza- 
tion amplitude is determined by the symmetry of the under- 



lying crystal-field Kramers doublets of the rare-earth ion and 
it is precisely this symmetry that is responsible for non-trivial 
topological structures in a KI. To analyze this topology, we 
first employ a tight-binding model on a simple cubic lattice, 
which is adiabatically connected to the Hamiltonian of the KI 
material. We show that in this case, band topology is uniquely 
determined by the noninteracting band structure of the system 
in the absence of hybridization. Secondly, we consider a more 
general KI in a lattice with a body centered cubic lattice, and 
show that regardless of microscopic details, there always ex- 
ists a parameter range in which a KI is a strong topological 
insulator (STI). 
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FIG. 1: Values of strong and weak topological indices and signs of 
Si (see text) at the high-symmetry points of the Brillouin zone (BZ) 
are shown as a function of position of the renormalized /-level rel- 
ative to the bottom of the conduction band. For the simple cubic 
tight-binding spectrum of the form 6k = —2t 2Za=x, y ,z cos ka topologi- 
cally strong, weak Kondo insulating behavior as well as regular band 
insulator (\sf\ > 6t) can be realized. These results are shown to be 
generic and independent of particular type of Bravais lattice. 

We begin with the periodic Anderson Kondo lattice Hamil- 
tonian, written in terms of the fermion operators associated 
with the crystal field symmetry of the underlying lattice 
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where £k is the dispersion of a tight-binding band of conduc- 
tion electrons. We assume that the ground state of the iso- 
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lated magnetic ion is a Kramers doublet |For), where T labels 
a particular representation of the crystal symmetry group. For 
instance, in a cerium-based Kondo system, the Ce ion has a 
Ce 3+ valence configuration, hence we have one /-electron in 
a J — 5/2 atomic shell. In Eq. (Q} above, the operator c \ a cre- 
ates an electron on site j in a Wannier state described by the 
quantum number a = ±, is the bare energy of the /-level, 
V is the bare hybridization, and Uf describes Hubbard repul- 
sion between /-electrons. One can relate the Wannier states at 
site j as follows JiSEitl: c ja = TjkA^>rk]aa-Cko-e' kRj , where 
the form factors [OpkW are two dimensional matrices 

[OrkW = (kTalkcr) = £ /r«|3m ,Lr\ ^(k). (2) 
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is a tight-binding generalization of the spherical Harmonics 
that preserves the translational symmetry of the hybridization, 
O(k) = <J)(k + G), where G is reciprocal lattice vector. Here, 
R are the positions of the Z nearest neighbor sites around the 
magnetic ion. 

The low-energy properties of the model (Q~|) are described 
in terms of renormalized quasiparticles formed via strong hy- 
bridization between the c- and /- states and on-site repul- 
sion Uf. In the regime where the / states are predominantly 
localized, we can neglect the momentum dependence of the 
/-electron self-energy £/(k, cj) - so that the effective 

low-energy Hamiltonian reads [ 12]: 
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where ^ is the bare spectrum of conduction electrons taken 
relative to the chemical potential, Sf = z[Sf + 2/(0)] is the 
renormalized /-level, V = \[zV , z = (1 — <9E/(w)/3<j)~L,, and 
1 denotes the unit 2x2 matrix. The KI is formed if the chem- 
ical potential of the quasiparticles lies inside the hybridiza- 
tion gap, separating the two bands with the spectra E ± (k) = 

\[& + s f ± A /(ft-e/) 2 + 4|VA k | 2 ],withA2 = iTr[^ k <D rk ]. 

From Eq. ©, we see that the form factors Opk are 
momentum-dependent unitary matrices that relate the spin 
quantization axes of the Bloch states and the spin-orbit cou- 
pled Wannier states. Our choice of hybridization ensures that 
the mean-field Hamiltonian (Eq. |4j is a periodic function sat- 
isfying U,/(k) = ^/(k + G). 

These form factors are uniquely determined by the wave- 
functions of a magnetic ion, \Ta). For example, in a tetragonal 
crystal field environment, the j = 5/2 Ce multiplet is split into 
three doublets: |rf ±) = | ± 1/2), |rf ±) = cosQ3)| + 3/2) + 
sinOS)|+5/2),and \Tf±) = sin(/3)|+3/2)-cos( J 6)| ±5/2), where 
the mixing angle j3 defines orientation of the corresponding 
states. In an orthorhombic environment, the Kramer's dou- 
blets are generally described by a linear superposition of all 



three wave-functions \Y {o) +) = u\ \ ± 1 /2) + v| + 3/2) + w\ ± 5/2) 
lfl3i [3. To discuss the topological properties of these KI 
we need to consider separately the form factors for dif- 
ferent T's. We distinguish these states according to their or- 
bital symmetry parameterized by the index a - 1,2,3 and the 
pseudo-spin quantum number (a = +) 11411 . Hence, we have 
//JO) = | ± 1/2), fl\0) = | ± 3/2), and ^J0) - | ± 5/2). 

The momentum-dependence of the hybridization gap 
A a (k) follows from Eq. (0. At small momenta k, 

Ai(k) = ^^[12008(20) + 5(3 + cos(40))] 1 / 2 , A 2 (k) = 

ljl\sm6\[n + 15cos(20)] 1 / 2 , and A 3 (k) = i^sin 2 6, 

where 9 and define the direction of the unit vector k, as- 
sociated with the point on the Fermi surface. Note that the 
hybridization gap has a line of nodes along the z-axis for the 
shapes a = 2,3, but generic combinations of all three form- 
factors characteristic of contain no nodes. The key results of 
this paper are most simply illustrated using the nodeless a = 1 
Kramers doublet as the ground-state of the magnetic ion. 

To analyze the topology of the bands we use the fact that 
topology is invariant under any adiabatic deformation of the 
Hamiltonian. We begin our study with a tight-binding model 
for a KI on a simple cubic lattice. The technical analysis is 
readily generalized to more complicated cases as discussed 
below. The most important element of the analysis is the odd 
parity form factor of the / electrons, <5 fl (k) = -O fl (-k). This 
parity property is the only essential input as far as the topo- 
logical structure is concerned. 

In Ref. 11511 . Fu and Kane demonstrate that in an insu- 
lator with time-reversal and space-inversion symmetry, the 
topological structure is determined by parity properties at the 
eight high-symmetry points, k* n , in the 3D BZ which are in- 
variant under time-reversal, up to a reciprocal lattice vec- 
tor: kJJ, = -k*, + G (see insets in Fig. Q}. In our case, 
these symmetries require that < ?4 ! /(k) = P'H m f(-k)P^ 1 and 
•7/ m/ (k) r = T'Hmfi-kyT' 1 , where the parity matrix P and 
the unitary part of the time-reversal operator T are given by 
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where o~2 is the second Pauli matrix. For any space-inversion- 
odd form factor, it follows immediately that <t> fl (k) = at a 
high-symmetry point. Hence, the Hamiltonian at this high 
symmetry point is simply < H mf (k* n ) = + s f )I/2 + (£ k ;„ - 
s/)P/2, where / is the four-dimensional identity matrix. 

The parity at a high symmetry point is thus determined by 
S m = sgn(4 r k;. - Sf). Four independent Z2 topological in- 
dices lHal (one strong and three weak indices) can be con- 
structed from 5 m : (i) The strong topological index is the 

8 

product of all eight 6 m 's; Tsti = Y\ S m = ±\; (ii) by set- 

m=l 

ting kj = (where j = x, y, and z), three high-symmetry 
planes, Pj = {k : kj = o}, are formed that contain four high- 
symmetry points each. The product of the parities at these 
four points defines the corresponding weak-topological index, 



3 



7 wti = ri dm - ±1 > j - x,y,z. The existence of the three 

k,„ePj 

weak topological indices in 3D is related to a Z2 topological 
index for 2D systems (a weak 3D TI is similar to a stack of 
2D Z2 topological insulators). Because there are three inde- 
pendent ways to stack 2D layers to form a 3D system, the 
number of independent weak topological indices is also three. 
A conventional band insulator has all of the four indices 

^STI = /wTI = 7 WTI = 7 WTI = +1 - While atl indeX 7 = 

indicates a Z2 topological state with the odd number of sur- 
face Dirac modes. For a KI with £k;„=o < e/ and > £/, 
we find /sti = ~~ 1 , and hence the Kondo insulating state is 
a strong-topological insulator which will be robust with re- 
spect to disorder. Weak-topological insulators and topologi- 
cal^ trivial insulators can in principle be found for different 
band structures and different values of s/, Fig. 1. A particu- 
larly interesting possibility is to tune topological phase transi- 
tions between different types of insulators (e.g., by applying a 
pressure), Fig. 1 . Although we have been specifically consid- 
ering a tight-binding model with a primitive unit cell, all our 
conclusions apply directly to systems adiabatically connected 
to this model. In order to prove it explicitly and to investigate 
more general cases, we develop a different and more general 
technique similar to that proposed in Ref. §71] . 

We shall now study an example of a KI for the specific 
shape a — 1 and explore the parameter range in which it re- 
mains a STL Here, the form factor is universally determined 
in the small-momentum limit by the /-wave symmetry of the 
electron orbitals and the point group symmetry of the lattice. 
Expression for the form factor of a Kramer's doublet F\ is 



O 



lk 
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where k = k/k and Y l m {V) are the spherical harmonics. At 
larger momenta, the form-factor depends on the microscopic 
details of the lattice and the Kondo coupling. In general form 
factor can be written as Oik = n\ ■ <x. 

We show below that momenta where Oik = are crucial 
in calculating the topological indices, so that our results are 
generic for a linear combination of all three shapes. The most 
obvious zero-point is located in the origin as shown in Eq. (O. 
We now prove that this zero-point is topologically protected 
and that its existence necessarily yields the existence of other 
zero-points Oik, = with k, + (this conclusion is sim- 
ilar to fermion doubling of relativistic fermions, which re- 
quires the presence of an even number of Dirac points on a 
compact manifold). To prove this, we draw a sphere S 2 at 
the origin with radius £0 as shown in Fig. |2b) and require 
Oik + on the sphere, which enables the definition of a 
Chern number C = ^ <j^, dS ■ e ,j *H/(Vk«j) x (Vk/ijO, where 

«; = mj Jm% + m 2 +m 2 and m,- is the component of 4*i 
as defined above. The Chern number is a topological index 
quantized to an integer value. For small ko, Eq. (O is asymp- 
totically accurate, which gives C = 1 . Notice that this topo- 
logical index is invariant as Icq changes adiabatically. The non- 



zero value of C indicates that £0 cannot be decreased to zero 
smoothly, and hence a point with Oik = must exist inside the 
sphere, which ensures Oik=o = 0. Since the BZ has a periodic 
structure and is a compact manifold, the same argument re- 
quires zero points Oik,„ = outside the sphere with |k|,„ > ko 
(m = 1,2...). This conclusion can be verified explicitly in the 
simple cubic lattice model studied above, where Oik = at 
each of the eight high symmetry points. 




FIG. 2: (a) Surface plot of the hybridization gap Ak for the form fac- 
tor of shape a = 1, Eq. (3). The ^-dependence of the hybridization 
gap originates from the underlying orbital structure of the localized 
electrons, which form the magnetic moment, (b) The BZ of an or- 
thorhombic lattice (or lattice with a higher point-group symmetry). 
The dots at the origin, face centers, bond centers, and the corner mark 
the eight high-symmetry points. The sphere at the origin separate re- 
gions I (inside) and II (outside). The intersection of the sphere with 
the high symmetry plane (k z = 0) is marked by the solid line. 

Now we relax the assumption about the simple cubic lat- 
tice and allow for a more general structure. Due to time- 
reversal and space-inversion symmetries, both bands are 
doubly-generated in a KI. Therefore, the corresponding Bloch 
wavefunctions *F'(k) and *F 2 (k) (which are four-component 
vectors) can be chosen arbitrarily up to a "local" U(2) trans- 
formation in momentum space. For fermions, T 2 = —I, and 
hence, we require that T^k) = ejj[ v P-'(— k)]* under time re- 
versal with eij being the Levi-Civita symbol. For concrete- 
ness, we focus below on the case with ^k < £/ at k — and 
£k, > Bf at all other zero points of Oik, which gave us a strong 
topological insulator in the model with a simple cubic lattice 
discussed above. For such a band structure, it can be easily 
checked that the wave-functions can not be defined globally 
in the entire BZ with the constraint T^k) = eypF(-k)]* 
lfl7ll . However, the BZ can be separated into two regions (c.f., 
Ref. dl) / (|k| < ko) and // (|k| > k ) with k being an ar- 
bitrarily small momentum such that only one zero point of 
Oik, k = is enclosed inside the sphere, Fig. |2jb). In each 
of the two regions, singularity-free Bloch wave-functions can 
be constructed. For example, the valence band has following 
wave-functions in region (/) 



¥j(k)=Al(0,£_,-Oi 2 ,-O22), 
^(k) = //_(£-,(), -On, -0 2i ), 

and in region (77) 

<(k) = // + (-a>2i.-^2.o.^). 
^ / (k) = jv + (-d>; 1 ,-*; 2 ,£ + ,o). 



(7) 



(8) 
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Here O^- are the (i, /(-components of the form-factor, 
£ ± (k) and Ak were defined above Eq. (|4j, and N± - 

l/^2V 2 A 2 k + (^k - Ef)E±. These two sets of Bloch wave- 
functions are connected by a "gauge" transformation *F)(k) = 
t/y(k)1'j / (k) at the boundary, S 2 , between the regions / and 
II, with the matrix U(k) = -a" z Oj . frjy. 

The topological structure of a 3D time-reversal invariant in- 
sulator is determined by the wavefunctions (|7|8t on the six 
high- symmetry planes Pj = Jk : kj - 0} and P'. = Jk : kj - 7rJ 
fiSIHtl- As shown in Fig.|2jb), the boundary between the 
two regions intersects with Pj and the intersections are cir- 
cles, C j. On such a circle, the matrix U takes the form of 
U (k) = exp{i(fo&Q + itpyi ■ &), with h being a fixed 3D unit vec- 
tor, corresponding to a specific "gauge." As a result, a winding 
number can be defined on C as Wj — § c * ■ Vk^k- However, it 
is defined modulo 2 only, because a gauge transformation can 
change tp^ — > ip^ + 2m<f> and hence change wj by an even num- 
ber 2m, where <p is the azimuth angle of k in a high symmetry 
plane (note that a transformation ip^ — > tp^ + (2m + \)<p would 
violate the symmetry constraint "FT" (k) = e, / [ v P 7 (-k)]*). For 
P'j, the corresponding winding numbers W, are zero for the 
case we studied here, because they do not intersect with the 
boundary. The topological indices can be computed from 
these winding numbers as follows lfl9ll : 



/STI = ("I)" 



and I J WTl 



(-D M 



(9) 



where 7sti is equivalently defined for j = x, y or z. For the 
types of KI band structures considered here, the topological 
indices can be universally determined by choosing a small 
enough ko and using the asymptotic form-factor of Eq. (O. 
As a result, we find that the generic Kondo system is a STI in 
full agreement with arguments above based on adiabatic de- 
formation of the Hamiltonian onto a simple cubic lattice. 

Let us briefly discuss the implications of our results for ex- 
isting Kondo insulators. From our theory's we expect that ma- 
terials in which /-electrons are close to integral valence are 
likely to be weak topological Kondo insulators and thus are 
unstable with respect to disorder. An interesting example is 
SmB6 for which recent LSDA+C/ band structure calculations 
112011 show the position of the /-level equals approximately one 
sixth of the bandwidth consistent with the core-level spec- 
troscopy measurements of the /-level occupation «/ ~ 0.7 
1 2111 . placing the quasiparticle f-level of SmB^ close to the 
border separating STI and WTI phases. Another promising 
candidate for the manifestation of topologically nontrivial in- 
sulating state is CeNiSn. Recent transport data in CeNiSn 
1 22[] shows suppression of semiconducting behavior in resis- 
tivity with increase in sample's quality, although there is an 
evidence for the gap formation at T 10K. Given that the 
/-electrons in these systems are predominantly localized, it is 
tempting to speculate that the CeNiSn is a weak topological 
Kondo insulator ascribing the semi-metallic transport proper- 
ties to metallic surface states. These are issues that we hope 



can be resolved in the near future through more accurate mod- 
elling and the use of high precision ARPES and STEM spec- 
troscopy. 

To summarize, we have developed a theory of topological 
3D Kondo insulators. Within our model topologically nontriv- 
ial insulating states are realized over a wide parameter range. 
In particular, we have shown that strong topological insulating 
state occurs when the position of the re-normalized /- level 
is near the top, or the bottom of the conduction band. This 
suggests the most likely candidates for this kind of behavior 
are heavy fermion materials which are more mixed valent or 
have narrow conduction bands. 
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